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£? Abstract 

CN We report a combined experimental and theoretical study of collision-induced dipolar relaxation 
> 

^^ in a cold spin-polarized gas of atomic nitrogen (N). We use buffer gas cooling to create trapped 

^T samples of ^'^N and ^^N atoms with densities (5 it 2) x 10^^ cm~^ and measure their magnetic re- 

C^^ laxation rates at milli-Kelvin temperatures. These measurements, together with rigorous quantum 

O 

^^ scattering calculations based on accurate ah initio interaction potentials for the ^S+ electronic 



cd 



state of N2 demonstrate that dipolar relaxation in N + N collisions occurs at a slow rate of ~10' 



13 



/\, cm^/s over a wide range of temperatures (1 mK to 1 K) and magnetic fields (10 mT to 2 T). The 

5-! 



calculated dipolar relaxation rates are insensitive to small variations of the interaction potential 
and to the magnitude of the spin-exchange interaction, enabling the accurate calibration of the 
measured N atom density. We find consistency between the calculated and experimentally deter- 
mined rates. Our results suggest that N atoms are promising candidates for future experiments on 
sympathetic cooling of molecules. 
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I. INTRODUCTION 

Owing to their unique controllability, cold and ultracold molecular gases hold promise 
for many research applications, ranging from quantum information processing [H [2] and 
simulation of condensed- matter systems |Ij to novel constituents of exotic quantum phases |3] 
and reagents for external field-controlled chemical reactions [4J. At present, cold molecular 
ensembles can be produced by a number of experimental techniques [H ISHTO]. which can 
be broadly classified as direct and indirect. The direct techniques are based on removing 
thermal energy from a pre-existing ensemble of molecules via collisional thermalization or 
time-dependent electromagnetic fields. Among the techniques of this kind are cryogenic 
buffer-gas cooling [S], velocity filtering [TU], and Stark and Zeeman deceleration [7HS]- The 
indirect cooling methods are based on creating molecules in ultracold atomic gases via 
photoassociation [HI [12] and sweeping a dc magnetic field across a Feshbach resonance [I] . 

While direct cooling methods produce molecules with initial temperatures between 50 
and 200 mK suitable for cold collision experiments [HI El |9], further cooling is required 
to reach the ultracold regime of interest to applications in condensed-matter physics and 
quantum information processing [1]. This may be accomplished via sympathetic cooling, 
a technique based on collisional equilibration of thermal energy, which takes place when a 
gas of molecules is brought into thermal contact with a cold reservoir of atoms. Because 
sympathetic cooling is driven by elastic collisions, it is a truly general technique, which has 
found numerous applications in cold atom physics [T314T5] . Most of the sympathetic cooling 
experiments performed so far used alkali-metal atoms (typically *^Rb [131 E]) because of 
their easy availability via laser cooling and their attractive collisional properties, which allow 
for sustainable evaporative cooling down to quantum degeneracy. 

The majority of experiments with cold molecules use permanent magnetic or electrostatic 
traps to capture molecules in their low-field-seeking Stark or Zeeman states, which are 
intrinsically unstable and may decay by collisions with background gas atoms. While elastic 
collisions lead to cooling, inelastic collisions heat the gas, cause trap loss, and shorten the 
lifetime of trapped molecules. The number of elastic collisions per inelastic collision must 
be large enough (>100) to allow for rapid thermalization while keeping inelastic losses to 
a minimum. An ideal atomic collision partner (X) for sympathetic cooling of magnetically 
trapped molecules must therefore meet the following acceptability criteria 



1. Be available in copious quantities at low and ultralow temperatures. In particular, 
atoms with magnetic moments of l/i^ or more (where fiB is the Bohr magneton) can 
be loaded in permanent magnetic traps via buffer-gas cooling [B] , evaporatively cooled 
to very low temperatures [I6], and co-trapped with molecular species [T7] . 

2. Have low inelastic X-X collision rates, so that sufficient density of X can be maintained 
in the trap at all temperatures to allow for efficient molecule thermalization; 

3. Have low inelastic collision rates with the diatomic molecules of interest, so that elas- 
tic atom-molecule collisions which drive thermalization occur more frequently than 
inelastic collisions. 

Previous theoretical work has addressed the optimal choice of atomic collision partners 
for sympathetic cooling of molecular species [T81 - I20] . The alkali- metal atoms, which satisfy 
the requirements (1) and (2), have so far received most attention and ab initio calculations 
of interaction energies and low-temperature collision properties have been reported for the 
alkali-metal atoms Rb and Cs interacting with OH [18] , NH |19] , and ND3 [20j . These studies 
have shown that the interactions between the alkali-metal atoms and diatomic molecules 
are strongly anisotropic, giving rise to large inelastic collision rates at low temperatures, 
thereby limiting the efficiency of sympathetic cooling [191 EO]- A recent theoretical study 
has shown that sympathetic cooling of OH(^H) molecules in low-field-seeking Zeeman states 
of e symmetry may be facilitated by superimposed electric and magnetic fields [21] . We note 
that certain trapping techniques employing ac electric [22], optical dipole [25] or microwave 
[23] fields, allow for trapping ground-state molecules, thereby eliminating the possibility of 
coUisional relaxation. At their present stage of development, however, these techniques are 
less advanced than magnetic or electrostatic trapping [1]. 

Recently, Wallis and Hutson theoretically explored the possibility of using the alkaline- 
earth atoms to sympathetically cool paramagnetic molecules [23]. Their ab initio calculations 
on Mg + NH collisions in the presence of an external magnetic field demonstrated that in- 
elastic collisions are suppressed at low magnetic field strengths, indicating that sympathetic 
cooling of NH molecules by collisions with laser-cooled Mg atoms might be possible [2^ . 
However, the detrimental inelastic collisions become very efficient at large magnetic fields, 
which requires precooling of molecular gas to very low temperatures (1 mK) before the sym- 
pathetic cooling can begin. In addition, the alkaline-earth atoms in their electronic ground 



states are not paramagnetic, which makes it challenging to produce large numbers of Mg 
atoms required for coUisional thermalization (criterion 1). The same limitation applies to 
recent proposals for using laser-cooled rare gas atoms to sympathetically cool large molecules 
like benzene [25] . 

We have recently suggested that molecular species bearing magnetic moments can be 
sympathetically cooled by collisions with spin-polarized nitrogen (N) atoms in a permanent 
magnetic trap [T71 [22] • Due to their large magnetic moments of S/i^, the N atoms can be 
efficiently confined in a magnetic trap using buffer-gas cooling [B]. Previous experimental 
work [17] has demonstrated that samples of N atoms with densities ~10^° - 10^^ cm~^ can be 
routinely produced and held in a magnetic trap for as long as 10 s, allowing for co-trapping 
of molecular species such as NH [17]. However, the cross sections for inelastic relaxation in 
N + N collisions were not measured in these preliminary experiments, leaving the question 
open of whether N + N collisional thermalization would prevail over two-body inelastic 
losses (criterion 2). While it is well-known that the two-body losses in doubly spin-polarized 
atomic gases are induced by the magnetic dipole interaction [27], the time scale for these 
processes in N + N collisions is unknown. In addition, the density of trapped N atoms could 
not be accurately determined due to the difficulties encountered in N atom detection. 

Here we present a combined experimental and theoretical study of low-temperature col- 
lisions in a cold spin-polarized gas of atomic nitrogen. We use buffer-gas cooling to trap 
large numbers of ^^N and ^^N atoms and study their collision-induced dipolar relaxation at 
milli-Kelvin temperatures. To interpret the experimental observations, we perform accu- 
rate ab initio calculations of the interaction potential between two spin-polarized N atoms 
and rigorous quantum scattering calculations of trap loss dynamics. By analyzing various 
sources of uncertainty in our theoretical results, we infer the upper and lower bounds to the 
calculated relaxation rates, which allows us to calibrate the N atom density based on the 
experimental measurements of trap lifetimes. 

The paper is organized as follows. In Sec. II, we describe our experimental apparatus 
and present measurements of collision-induced trap loss rates for both ^'^N and ^^N isotopes 
of atomic nitrogen. Sections IIIA and IIIB present ab initio calculations of the interaction 
potentials for N2 and give a brief outline of quantum scattering calculations on N + N colli- 
sions. Sections IIC and IIIC compare our theoretical results with experimental data. Section 
IV gives a brief summary of our results and outlines possible future research directions. 



II. EXPERIMENT 



A. Apparatus 



The experimental apparatus is similar to those described in Refs. [T71 [2S]. A diagram 
of the trapping region is shown in Fig. [ij The trapping region is centered about a pair 
of super-conducting magnetic solenoids that produce a spherical-quadrupole magnetic trap 
with depth of 3.8 T. For atomic nitrogen, with a magnetic moment of 3/iB, the corresponding 
trap depth is about 7.6 K. The magnetic field contours are shown in gray in Fig. [TJ In the bore 
of the super-conducting solenoids resides a cylindrical copper buffer gas cell, maintained at a 
temperature of about 600 mK by thermal anchoring to a ^He refrigerator. The buffer gas cell 
has an aperture at both ends of the trapping region. At one end, a 1 cm diameter aperture 
allows the atomic nitrogen to enter the trapping region from a room temperature atomic 
beam. At the opposite end, a 3.80 cm diameter aperture allows the buffer gas to be rapidly 
introduced into and subsequently removed from the trapping region. The atomic nitrogen 
beam is generated using a DC glow discharge source with N2 as the process gas, operating at 
a stagnation pressure of 100 torr . The atomic source is turned on for approximately 40 ms 
to load atoms into the trapping cell. Simultaneous with the introduction of atoms into the 
trapping region, '^He buffer gas is introduced into the cell by pulsing open the cryogenic 
buffer gas reservoir [28j. The density of the buffer gas during loading of the atoms is on the 
order of 10^^ cm~'^. The loading pulse of buffer gas then exits the buffer gas cell via the 
3.80 cm aperture on a time scale of 50 ms. Following loading of N into the magnetic trap, the 
final background buffer gas density in the trapping region is approximately 10^^ cm~^. The 
density is set by the rate of helium desorbing from the buffer gas cell walls [29] . Nitrogen- 
helium collisions, occuring at a rate of about 100 Hz due to residual buffer gas, pin the 
trapped nitrogen temperature to the temperature of the cell walls. This allows the trapped 
nitrogen temperature to be monitored using a ruthenium oxide thermometer mounted to 
the cell wall. A resistive heater mounted to the cell wall allows the nitrogen temperature 
to be adjusted. The nitrogen trap loss lifetime due to these N-He elastic collisions is on the 
order of 100 s since the temperature of the He buffer gas is more than a factor of 10 lower 
than the depth of the magnetic trap. This long trap lifetime makes it possible to study N-N 
collisions that lead to trap loss on time scales of 10 s. 



B. Atomic nitrogen detection 

To detect the trapped atomic nitrogen we use two-photon absorption laser induced fluo- 
rescence (TALIF) [30]. We excite atomic nitrogen in the ground (2^^)^83/2 state by absorbing 
two photons at 206.7 nm to the excited {3p)^S3/2 state at 96750 cm~^. The excited (3p)''S3/2 
state has a hfetime of 26 ns [31j and decays to the (3s)^P states, emitting hght near 745 nm. 
Aim focal length lens placed outside the vacuum chamber is used to focus the excitation 
laser onto the trapped sample. The fluorescence is collected using a lens mounted at the 
midplane of the magnet and sent to a photomultiplier tube for detection. 

Estimation of the trapped atomic nitrogen density from the TALIF signal is difficult. 
Both the fluorescence collection efficiency and nitrogen excitation probability are required 
to convert the TALIF signal to an absolute nitrogen density. From geometric considerations 
and fluorescence measurements using trapped NH, we estimate our fluorescence collection 
efficiency to be 10~^ [26]. This value for fluorescence collection effieciency we estimate to 
be accurate to within a factor of 3. To calculate the nitrogen excitation probability, one 
needs precise knowledge of the spatial, temporal, and spectral properties of the excitation 
laser. In the low intensity limit, where depletion of the ground state and photo- ionization 
are negligible, one can show that the total number of fluorescence photons produced from 
the sample is [32j: 

iVphoton = ^^'^7^ / n^r{r)S\r)dV / F\t)dt (1) 

where A^photon is the total number of fluorescence photons produced per pulse, a^'^^ is the 
effective two-photon cross section , ngr(^) is the nitrogen ground state density distribution, 
E is the laser pulse energy, hu is the excitation laser photon energy, and S{r) and F{t) 
are the normalized spatial and temporal proflles of the laser beam {J S{r)dA = 1, and 
J F{t)dt = 1). Here we have assumed the spatial and temporal variations are independent. 
We can express the value of the spatial integral in terms of an effective 1/e^ beam waist Wq, 
where J S{rYdA = ti^^Wq^. Similarly, the temporal integral can be expressed in terms of an 
effective pulse duration, Tgx, where f F{tYdt = a/2 ln(2)7r^^(rcx)^^. The resonant effective 
two-photon cross section can be expressed as [33] : 



a 



(2) = ^(2)^(5 = 0)G(^)(i = 0) (2) 



where g{5 = 0) is the resonant hne shape factor, G^'^\t = 0) is the second-order intensity- 
correlation function of the excitation laser, 6 is the detuning from the atomic resonance, 
6 = 2u;iascr — ^0, and the line shape factor is normalized such that j g{6)d6 = 1. 

The spatial profile of the laser at the position of the atoms is measured using a CCD 
camera. The value of J S{rydA is measured to be 7r~^(120/im)~^. The term J ngT.{r)S'^{r)dV 
in Eq. ([I| can be evaluated in the following manner. For the trap geometry in our experiment, 
the trapped nitrogen density does not vary significantly over the spatial profile of the laser. 
In the direction of propagation of the laser, the nitrogen atoms are confined to a effective 
length, les = 2 mm. The term can then be evaluated J ng^{r)S^{r)dV = riolcs J S{r)'^dA, 
where hq is the nitrogen density at the center of the trap. 

To monitor the temporal profile of the laser we pick off a portion of the laser beam and 
direct it onto a ceramic beam dump. We then use a fast photodiode to monitor the light scat- 
tered from the beam dump. The value of J F{tydt is measured to be A/21n(2)7r~i(9.5ns)~^. 
The photodiode is also used to monitor the laser pulse energy, E, by cahbrating its signal 
using a pyroelectric energy meter. 

The spectral profile of the pulsed laser is more difficult to characterize. The commercial 
Sirah pulsed dye laser uses a Littman-Metcalf configuration for the resonator cavity [31] 
to produce laser light at 620 nm. The resonator cavity has a linewidth of about 1.5 GHz 
and longitudinal mode spacing of about 600 MHz. For each laser pulse, several different 
longitudinal modes may lase. This is observed in our setup using a Fizeau interferometer 
[35| |36] as a spectrum analyzer. The 620 nm light is subsequently frequency doubled in 
a KDP crystal, and the doubled light mixed with the fundamental in a BBO crystal to 
produce light at 207 nm. The shot to shot spectral variation of 1.5 GHz at 620 nm leads 
to a variation of 9 GHz (0.3 cm~^) at the resonant atomic frequency. For comparison, the 
expected Doppler broadening (full-width half-max) of the atomic transition is expected to 
be 300 MHz at 600 mK with Zeeman broadening of less than 100 MHz [371138] . As a result, 
the value of g{0) will be determined by the spectral properties of the laser. An upper limit 
on g{0) can be determined by directly measuring the observed nitrogen signal linewidth. 
Figure |2] shows a trapped nitrogen spectrum taken at a cell temperature of 600 mK. To 
acquire these spectra, we monitor the wavelength of the excitation laser using a wavemeter 
with resolution 0.01 cm~^ at 620 nm, corresponding to a resolution of 0.06 cm~^ at the atomic 
transition frequency. Each data point represents the average signal of three consecutive laser 



excitation pulses. These consecutive laser pulses have a corresponding frequency jitter on the 
order of 0.3 cm~^ at the atomic transition frequency. The measured linewidth of 0.76 cm~^ 
(23 GHz) is a result of a combination of the actual spectral linewidth of the laser and the 
limited resolution of our measurement technique. For calculations in this paper, we take 
^(0) = (2/7r)(27r x 10 GHz)-^ This value of ^(0) should be good to a factor of 2. 

Finally, due to the complicated spectral mode structure of the light at 620 nm, it is 
not clear what the value of G^'^\t = 0) is for the light produced at 207 nm. In Ref. |33] . 
Bamford and coworkers analyze values of G''^^^(t = 0) for various pulsed laser systems. They 
find G'^^^(O) ranges between 1.5 and 3.0 for typical multimode pulsed laser systems, though 
our setup is not directly measured in the reference. Taking a value of G^^\0) = 2 will be 
within 30% of the actual value. 

The parameters to calibrate the TALIF signal are summarized in Table I. Due to the large 
uncertainties in the spectral properties of the excitation light and the fluorescence collection 
efficiency, the nitrogen densities calculated using the values listed in Table I should be 
accurate only to within an order of magnitude. Also omitted from this analysis is the role 
of laser polarization and atom orientation, which would likely introduce corrections of order 
unity to the calculation. 

C. Experimental results 

Here we present our observations of trapping of atomic nitrogen and discuss the nature 
of the observed trap loss. In particular, we are interested in the ratio of the elastic N + N 
collision rate to the inelastic N + N collision rate, 7. Measurements of both the elastic and 
inelastic collision rates are desirable, though to measure each independently, one needs to 
have an absolute calibration of the atomic density. Since our estimates of atomic N density 
from the TALIF signal are only good to about an order of magnitude, we lack the information 
we need to make precise measurements of the elastic and inelastic N + N collision rates. 

However, it is possible to directly measure 7 without precise knowledge of our atomic 
nitrogen density by investigation of the dynamics of the magnetic trap loss. Qualitatively, 
for very deep traps, evaporation of the sample is suppressed, and loss is driven by inelastic 
N + N collisions. For lower trap depths, evaporative loss due to elastic N + N collisions 
can dominate trap loss. By measuring trap loss due to N + N collisions over a range of 
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trap depths, it is possible to directly extract 7. A discussion of our magnetic trap dynamics 
follows. 

The expression for loss of atomic N in our magnetic trap has the form (see Appendix A) 

no = -—-{Ki^^Ml no (3) 

7.6 THe 

where no is the peak trap density, (-fCioss) is the trap average 2-body loss rate coefficient, 
1/7.6 is a factor for our trap geometry, and thg is the 1/e lifetime associated with loss due 
to atom-helium collisions, (i^ioss) includes trap loss both from atom evaporation, {^fi'evap), 
and loss from atom-atom inelastic collisions, (-ft'm)- The evaporative portion of trap loss 
can be expressed as (-ft'evap) = f{v){^ci), where /(r^) is the fraction of elastic collisions that 
lead to atom loss at trap depth t] = Ttrap/^atom, where Ttrap is the trap depth expressed in 
units of temperature, and Tatom is the temperature of the trapped atomic sample. For our 
trap geometry, Monte-Carlo simulations of trap dynamics yield /(r/) = 1.9(?7 — 3) exp(— r/), 
which agrees well with analytic expressions for /{rj) [39j. The relationship between (-f^in) 
and the dipolar relaxation rates derived from the quantum scattering calculations in Sec. 



Ill is given in Appendix A. Combining these expressions, we have: 



Kloss) = fiv){Kol) + {Kin) (4) 

= (Kei>(/(r/) + l/7) (5) 

Equation ([s]) provides an expression for extracting 7 without precise knowledge of the ab- 
solute atomic nitrogen density. 

Figure [3] shows a typical nitrogen trap decay. For each data point in Fig. |3] we load N 
atoms into the trap at t = and detect the remaining N atoms after waiting a period of time 
between 2 s and 100 s. Attempts to continually detect trapped N atoms during a single trap 
loading result in rapid N trap loss, most likely due to optical pumping. The shot to shot 
variation in N signal is likely due to the variation of the spectral properties of the excitation 
laser. We fit our nitrogen time decay data to the solution of Eq. ^ to arrive at values for 

(i^loss)- 

Values for (ii'ioss) are measured for values of i] between 10 and 14, which correspond to 
Tatom between 550 niK and 650 mK with magnetic trap depths between 3.3 T and 3.9 T. A 
plot of (-ft'ioss) versus 77 is shown in Fig. 4 The solid line in Fig. 4 is the fit to Eq. fc| with 
(/^ei) and 7 allowed to vary as fit parameters. This yields a value of 7exp = (6.5 ± 5.5) x 10^. 



Figure |5] shows the results of quantum scattering calculations described in the following 
sections, which give 7theory = (1-0 ± 0.3) x 10^ at T = 600 mK. The dashed line in Fig. E^ 
shows a fit of Eq. ^ with the value of 7 fixed to 7theory = 1000, where the only free fitting 
parameter is (i^ci)- The values of 7exp and 7theory are consistent, and the value of 7 > 1000 
is favorable for evaporative cooling of atomic N in a magnetic trap. 

From the data in Fig. [4], one can calibrate the trapped atomic nitrogen density by fitting 
the observed loss rates to the calculated loss rates. What is needed is simply a scaling factor, 
uq = CN^obs, to convert the observed nitrogen signal, riobs, to an actual nitrogen density, no- 
Substitution of the scaling factor relation into Eq. [3] results in the expression: 

^obs = -;r77(-ft'loss)cN'^obs '^obs (6) 

7.6 THe 

By setting the value of (/Cioss) to the theoretical value calculated as described in Sec. Ill 
and allowing cn to vary as a fit parameter, one can then fit the solution of Eq. [6] to the 
observed nitrogen trap decay to arrive at a value for the scaling factor cn- When using 
this calibration technique we need to consider the systematic uncertainties associated with 
arriving at theoretical value of (-ft'ioss)- In particular, for 77 > 12, (-ft'in) accounts for 60% or 
90% of the total (-ft'ioss) for values of 7 = 6500 or 1000, respectively. Although there may 
be large uncertainty in the actual value of 7, for large r], the uncertainty in the systematic 
correction of (-ft'm) to arrive at a total (-ft'ioss) is only about 30%. From this calibration 
method we estimate we typically trap atomic nitrogen at initial peak densities of more than 
(5 ±2) X 10^^ cm~^, corresponding to more than (3± 1) x 10^^ trapped nitrogen atoms. The 
error in these nitrogen density measurements is dominated by the systematic uncertainties 
associated with the model of trap loss dynamics and the quantum scattering calculations 
of (^Kioss), but also includes the statistical uncertainties associated with the experimental 
measurement of (-ft'ioss) (or equivalently cn) . The atomic nitrogen density calculated using 
this method agrees with our nitrogen density estimates from the TALIF signal to within 
a factor of 5, consistent with our expected uncertainty in the TALIF signal of an order 
of magnitude. This technique for measurement of the trap nitrogen density is valuable 
due to technical difficulties associated with a direct spectroscopic measurements of trapped 
atomic N. 

We also observe trapping of the bosonic isotope ^^N by using isotopically enriched (> 
98%+) ^^N2 as the process gas. No differences in trap loss were observed between ^^N and 
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^^N at a trap temperature of 600 mK. 

III. THEORY 

A. Ab initio calculations of interaction potentials 

To evaluate the potential energy curve (PEC) for the '^S^ electronic state of N2, we use 
two different ab initio approaches. The first approach is based on the coupled cluster method 
including perturbative triple excitations (CCSD(T)) and extrapolation to the complete basis 
set (CBS) limit. The second approach uses the coupled cluster method with full iterative 
triple excitations (CCSDT) and a fixed basis set with additional bond functions. The inter- 
action energies in both methods are calculated within the supermolecular approach, where 
dimer and monomer energies are calculated with dimer centered basis sets. We applied Boys 
and Bernardi counterpoise procedure to correct for the basis set superposition error. 

In the first approach, the PEC was calculated using a single-reference restricted Hartree- 
Fock (RHF) wave function as a starting point, followed by a spin-unrestricted coupled cluster 
treatment jl3] with single, double and non-iterative triple excitations (UCCSD(T)) as im- 
plemented in the MOLPRO suite of programs [44] . The use of the single-reference approach 
is justified because the high-spin electronic state '^S+ can be well described by a single deter- 
minant. For the purpose of extrapolation to the CBS limit, we used a series of augmented, 
correlation-consistent triple-, quadruple-, quintuple- and sextuple- zeta basis sets of Dunning 
et al. HHIig denoted as AVTZ, AVQZ, AV5Z and AV6Z, respectively. The Is orbitals of N 
were frozen in these calculations. The interaction energies were calculated at 50 internuclear 
separations from R = 2.5 ao to -R = 50 qq and fit to analytic functions of R using the repro- 
ducing kernel Hilbert space (RKHS) method j^HTj. The radial kernel was composed of the 
short range part and the asymptotic long-range part proportional to R^^. The smoothness 
of the one-dimensional kernel parameter was set to 2 to allow for a smooth extrapolation of 
ab initio data points to the asymptotic region {R > 50ao)- 

For each R point we performed extrapolation to the CBS limit using AVTZ, AVQZ, 
AV5Z and AV6Z interaction energies. To fit the series of interaction energies we used the 
empirical formula Ex = Ecbs + Ae"'-^^^^ + Be^^^^^^ suggested by Peterson and coworkers 
[15| HH], where X = 3,4,5,6 is the number of "zetas" in the basis set. The resulting 
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UCCSD(T)/CBS PEC (labeled as potential A) is shown in Fig. pi The potential A has a 
minimum at R^ = 7.21 a^ with a well depth of Df. = 29.3 cm~^. 

The "^S^ state has a large multiplicity, leading one to expect a significant contribution 
of higher excitations in the CCSD(T) method. In order to estimate this contribution, we 
included the full iterative triple excitations in our ab initio calculations of the interaction 
energy for the ^S^ state. The inclusion of full connected triple excitations makes the ab initio 
calculations much more computationally demanding, and we employed a single correlation- 
consistent AVTZ basis set with an additional set of 3s3p2d2flg bond functions (BF) placed 
at the middle of the N2 bond to reduce computational costs. To perform the full MRCCSDT 
calculations, we used the MRCC program |19] by Kallay et al. [50] interfaced with the 
MOLPRO code P]. The MRCCSDT calculation used a single-determinant RHF wave 
function as a reference. The resulting MRCCSDT/ AVTZ+BF PEC (labeled as "Potential 
B") is shown in Fig. ^ The minimum of the potential B has a well depth of D^, = 31.6 cm~^ 
and is located at Re = 7.18 oq. These values may be compared with the previous ab initio 
results Re = 7.5 oq and De = 21 cm^^ obtained by Partridge et al. [51]. We note that with 
the same AVTZ+BF basis the calculated well depth at the UCCSD(T) level is similar to 
that of potential A. The inclusion of the full triple excitations thus increases the well depth 
by approximately 7%. 

To compare the asymptotic behavior of potentials A and B, we fitted their long-range 
parts between R = 20ao and R = 40ao to the form —Cq/R^ — Cs/R^. The fitting procedure 
yields Cq = 23.36 E^al for potential A and Cq = 24.0 E^al for potential B. These results 
are in close agreement with the highly accurate value of 24.2 [52], thereby attesting to the 
accuracy of our ab initio calculations. 

Table II hsts the bound levels of ^^N2(^E+) and ^^N2(^S+) calculated using potentials A 
and B in the absence of a magnetic field. Both potentials are deep enough to support three 
bound levels with f = 0, 1, and 2. The number of rotational levels decreases from 10 (or 11 
for potential B) for f = to two for v = 2. We note the presence of accidental degeneracies 
between the rotational levels corresponding to different v. The levels v = 0,i = 11 and 
V = l,i = Q calculated with potential A have very similar binding energies of —0.221 cm~^. 
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B. Scattering calculations 

The Hamiltonian for two S-state atoms such as N('^ 5*3/2) colhding in the presence of a 
uniform magnetic field of strength B may be written {h = 1) 

" = -2iiW|j« + 2U ^ "'"<"> ^ ^'"^"^ ^ "' ^ "^' <'* 

where i/ = A , B enumerates the atoms, /i is the reduced mass of the N2 molecule, £ is the 
orbital angular momentum for the collision, and R = \R\ is the interatomic separation. The 
Hamiltonian of the isolated atom u is given by 

H. = iX ■ S, + 2^l^BS,^ - y^4 , (8) 

where 5",^ and I^ are the electron and nuclear spins, /io is the Bohr magneton, 7j, is the 
hyperfine constant, and ^jy is the nuclear magnetic moment. In this work, we consider 
both naturally occurring nitrogen isotopes, fermionic ^^N [I^, = 1/2, 7/27r = 10.451 MHz) 
and bosonic ^^N [I^, = 1, 7/27r = 14.646 MHz) [53]. The operators S,^^ and ly^ yield the 
projections of 5*1, and I^, on the space- fixed quantization axis defined by the external magnetic 
field. The magnetic dipole-dipole interaction is 



where a is the fine-structure constant, [Sa ® 5'b] is a second-rank tensor product of atomic 
spin operators, and Ykq{R) are the spherical harmonics. The vector R = R/R describes the 
orientation of the N2 collision pair in the space-fixed coordinate frame. 

The spin-dependent interaction potential between the atoms may be written as 



V^''{R)= Yl E VsmSMs){SMs\, (10) 

S=\Sa-Sb\Ms=-S 

where 5* = 5'a + ^b is the total spin of the collision complex and Ms = Msj^ + M^^ is the 
projection of S on the space-fixed quantization axis. In this work, we use the accurate ah 
initio interaction potentials for the '^S^ electronic state of N2 {S = 3) calculated as described 
in Sec. IIIA. 



Equation (10) is parametrized by four spin-dependent interaction potentials of the N2 
molecule correlating with the lowest dissociation limit N^^ 6*3/2) + N(''S'3/2). In addition to 
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the high-spin '''S+ potential described above, the low-spin electronic states of ^S+, ^Sj, and 
^S^ symmetries (5* = 2, 1, and 0) should be taken into account. While the X^T,^ and 
A'^Tig electronic states were subject to several theoretical studies [S^IES], no high-quality 
ab initio calculations are available for the A'^'Et electronic state. In the absence of more 



accurate information, we choose to parametrize the spin-dependent interaction (10) by the 
Heisenberg Hamiltonian [571459] 



V^'^(R) = V^'iR) - 2JiR)SA ■ 5b, (11) 

where V^^ is a spin-independent interaction potential, and J{R) is the spin-exchange (SE) 



coupling strength. It follows from Eq. ( 11 ) that the interaction potentials for the spin states 



5* and S — 1 differ exactly by twice the SE coupling strength 

Vs{R)-Vs-i{R) = 2SJ{R). (12) 



Equation (12) allows us to obtain the four PECs in Eq. (10) in terms of two parameters: 
(i) the potential energy curve for the "^S^ state calculated in Sec. IIIA and (ii) the SE 
coupling strength J{R). For the latter, we use the expression derived by Smirnov and 
Chibisov [57HS9] 

J{R) = Ci?°e"^^, (13) 

where (3 = ySI and a = 7/(3 — 1 are expressed via the atomic ionization energy I (0.53412 
Efi for the N atom): (3 = 2.0671 a^^ and a = 2.3864. In Sec. IIIC, we will use C as a free 
parameter to vary the magnitude of the SE coupling in order to explore the sensitivity of 
scattering cross sections to the interaction potential. 

If the weak hyperfine interaction of ^^N is neglected (see Appendix B for a justification), 
the wave function of the N2 collision complex can be expanded in direct products of electronic 
spin functions and partial wave states 

where 



kV 



\SAMs^)\SnMs,)+v{-Y\SBMs,)\SAMsMime), (15) 



^Ms,Ms,em, 12(1 + 6ms,MsJV/'' 

IS^Ms^) are the electronic spin basis functions of individual atoms A and B, Ms^ are the pro- 
jections of Si, on the space-fixed quantization axis, and \imi) = YimdR)- The direct-product 
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basis (15) is properly ordered (ta > tb) and symmetrized to account for the quantum statis- 
tics of indistinguishable bosons (^^N, rj = 1, odd i) or fermions (^^N, rj = —1, even i). 



The matrix elements of the spin-dependent interaction potential (10) in the symmetrized 
basis are 



w 



sdi i»? 



Sff'S. 



rriim^ 



[(1+5- 



Ms^MsJi^ + ^M's^M'^^] 



11/2 



sdi 



X [(5aM5J(5bM5JV^^15aM' )|5bM' ) + r/(-)^(5AM5j(5BM5j\/^^|^BM' )|5aM' )]. 



(16) 



The second term on the right-hand side (which arises from symmetrization) can be obtained 
from the first term by interchanging the indices Mg^ -H- M'g^. We therefore only need to 
evaluate the unsymmetrized matrix element 



^Msj,Ms^emr,M's^M's^e'm'^ 



Sii'Sm,m',{SAMs^ I {SbMs,\V'''\SaM'sJ\SbM'sJ. (17) 



where the subscripts Sa, Sb and 5*^, S'^ have been omitted for clarity. Expanding the 
product of two spin functions in a Clebsh-Gordan series [UCTl lUT] , we obtain 



V, 



sd 
Ms, Ms„£m^;Af' M'i'm' 



Su'Sm.m', ^(25+l)(-)^^« 

S,Ms 



Sa Sb 
Ms^ Ms, 

Sa Sb 




Vs{R), (18) 



Because the spin-dependent interaction potential (10) is diagonal in the total spin S and 



its space-fixed projection Ms, the matrix elements between the fully spin-polarized initial 
state \Sa,Msj^ = Sa)\Sb,Msj^ = Sb) and all other spin states vanish identically. Thus, in 
the absence of the magnetic dipole interaction (see below), spin-polarized atoms can only 
undergo elastic scattering. 

The matrix elements of the magnetic dipole interaction can be derived as described else- 
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where [Ml ESI EH- Here, we only present the final result 



T/dip _ V30a ^ .jne+SA+SB-MsWop I TWO/ I i\ll/2 

g 

X [(25a + 1)5a(5a + 1)]^/'[(25b + 1)5b(5b + 1)]^/^ | ^ ^ ^' 

\^0 

^ ^ / 1 1 2 W 5a 1 5a 




Unlike the interaction potential (10), the magnetic dipole interaction does not conserve 



Ms and couples the fully spin-stretched state |5a,M5'^ = 5A)|5B,M5'g = 5b) to other spin 



states, thereby inducing spin-fiipping transitions. By transforming Eq. (19) to the total spin 
representation, one can show that the matrix elements of the magnetic dipole interaction 
vanish identically unless 5 — 5' = 0, ±2 [631 169]. Thus, the mechanism of dipolar relaxation 
in collisions of fully spin-polarized N atoms is completely determined by the electronic states 
of '^T,^ and '^S^ symmetries [69]. This mechanism is typical of light 5-state atoms [62]. The 
rate constants for dipolar relaxation can be measured by observing collision-induced loss of 
atoms from a magnetic trap as described in Sec. II. 

A system of close-coupled (CC) equations for the radial functions Fm^ Ms imX^) results 
when the expansion (14) is combined with the Schrodinger equation with Hamiltonian ([T]). 



The CC equations are integrated numerically on a radial grid extending from R = 4.0 ao to 
R = 50.0 Oq with a constant spacing of 0.04 Oq using the improved log-derivative algorithm 
|68j . The calculations are carried out separately for each total angular momentum projection 
M = Msj^ + Msj^ + Tn^, which is conserved for collisions in external fields. The scattering 



5-matrix is computed directly in the uncoupled basis (15) and used to evaluate the cross 



sections for collision-induced transitions between different Zeeman states [651 - 167] 

7r(l + OMs.Ms ) v-^ v-^ v-^ 

^Ms^Ms^^M's^M's^ = p —l^l^l^ \SMs^M'sjMs^M'sJie'Sm,m'^ 

^S^Msj^ M e,mee',m'i, 

~ ^Msj^Ms^tmf,M's M'g e'm'J ) (20) 

where the factor (1 + Sms Ms ) accounts for the indistinguishability of colliding atoms. In 
order to make sure that our numerical results are correct, we repeated scattering calculations 
with a different code [69] and obtained the same results. A total of 64 channels were coupled 
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for M = 0, corresponding to the partial wave states with £ < £max = 13. The calculated 
elastic and inelastic cross sections were converged to better than 3% with respect to the 
initial and final integration distances and imax for all collision energies and magnetic fields 
studied in this work. 



By averaging the symmetrized cross sections (20) over a Maxwell-Boltzmann distribution 
of collision energies in the interval Ec = 0.01 — 2 cm~^ with a grid step size of 0.01 cm~^, 
we obtain state-resolved dipolar relaxation rates as functions of temperature T 

Kms^Ms^^m'.^m',^ (T) = (^) I PEcaMs,Ms,^M',^M',^ iEc)e-^o/ksT^E^^ (21) 
where /5 = l/ksT and ks is Boltzmann's constant. 

C. Theoretical results 

1. Cross sections and rate constants for dipolar relaxation 

Figure [t] shows the energy levels of ^^N as functions of the applied magnetic field. At zero 
field, the ground state is split by the hyperfine interaction into three levels with F = 5/2, 
3/2, and 1/2. The inset of Fig. [T] shows the hyperfine splittings as functions of the applied 
magnetic field. As shown in Appendix A, the hyperfine splittings have a minor effect on 
N + N collisions except at very small magnetic fields (<20 G), so it is a good approximation 
to consider bare spin states \Sa,Ms^) with Ms^ = —3/2, . . . 3/2. The magnetic trapping 
experiments described in Sec. II select N atoms in the fully spin-polarized state IS* = Ms = 
3/2), so in the following we will only consider collisions of N atoms initially in this state. 

The calculated cross sections for elastic scattering and dipolar relaxation in ^^N + ^'^N 
and ^^N + ^^N collisions are displayed in Fig. |8]as functions of collision energy Ec- In the 
limit of vanishing collision energy, the variation of the cross sections with Ec is determined 
by the Wigner threshold law; the cross sections for elastic scattering vary as E"^ and those 
for dipolar relaxation as E(7 . Collisions of identical fermions such as ^^N are determined 
by p-wave scattering, so both the elastic and inelastic cross sections vanish as Ec tends to 
zero. The situation for the bosonic isotope ^^N is the opposite: the inelastic cross section 
diverges and the elastic cross section approaches a constant. At collision energies above 
~0.01 cm~^, partial waves with i > 1 begin to factor in, altering the dependence of the cross 
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sections of collision energy, and leading to the appearance of shape resonances. A particularly 
pronounced resonance at Eq ~ 1.3 cm"^ corresponds to the quasibound ^^N2(^S+) molecule 
in the C. = 7 rotational state. 

In Fig. [9} we plot the cross sections for dipolar relaxation as functions of collision energy 
and magnetic field. Magnetic fields increase the splitting between the incoming and outgoing 
collision channels and suppress dipolar relaxation [63] . The suppression is quite pronounced 
at high magnetic fields (on order of 1 T), but does not alter the dependence of the cross 
sections on collision energy. In particular, increasing the magnetic field from 0.01 T to 1 T 
reduces the lifetime of the £ = 7 shape resonance by more than two orders of magnitude, 
but leaves its position intact. The same is true for the low-energy £ = 2 shape resonance in 
^^N. As shown in Fig. [9] for both ^^N and ^^N, the decrease of the cross sections with B is 
not always monotonous. 

As discussed in the Introduction, collisional stability is a key ingredient to efficient use 
of trapped atomic gases for sympathetic cooling of molecular ensembles. Fig. [Io|^a) shows 
the rate constants for dipolar relaxation calculated by thermally averaging the cross sections 
shown in Fig. [8] at a fixed magnetic field of 0.1 T. The relaxation rates for both N isotopes 
display broad maxima at T ~ 10 mK (for ^^N) and T ~ 50 mK (for ^^N). The rate for 
^^N features an additional maximum near 1 K due to the £ = 7 shape resonance shown in 
Fig. |8| At temperatures below 5 mK, inelastic collisions occur in the Wigner s-wave regime, 
and the rate constants for dipolar relaxation tend to zero for ^^N and approach a constant 
value of 5.5 x 10^^'^ cm^/s for ^^N. The ratios of the rate constants for elastic scattering and 
dipolar relaxation displayed in Fig. [lO[b) remain high (7 > 100) down to ~10 mK for ^^N 
and ~2 niK for ^^N. This result shows that trapped ensembles of ^^N and ^^N atoms with 
densities 10^^ cm~^ will have lifetimes ~2 s over a wide range of temperatures from 1 mK 
to 1 K. We note that because the elastic cross section for ^^N becomes very small at T < 1 
mK, ^^N is a more promising candidate for cooling molecules to temperatures below 1 mK, 
whereas both ^^N and ^^N isotopes appear suitable for sympathetic cooling to temperatures 
above 1 mK. 



Figure 11 shows the temperature dependence of state-to-state rate constants for dipolar 



relaxation in ^'^N + ^^N collisions (21). The single spin-fiip transition jM^^ = 3/2)|M5g = 
3/2) —7- \M'g^ = l/2)\Mg^ = 3/2) dominates over the whole range of magnetic fields at 
both 0.1 K and 0.6 K, and the double spin-fiip transition \Msj, = ?>/2)\Ms^ = 3/2) -^ 
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\Mg^ = l/2)\Mg^ = 1/2) is the next most efficient. The rate constants for other transitions 
(changing Ms by 2 or more) are neghgibly small. 

2. Sensitivity of trap loss rates to the interaction potential 

In order to verify the reliability of our scattering calculations, it is essential to analyze 
various sources of uncertainty that can affect the accuracy of numerical results for the 
dipolar relaxation cross sections and trap loss rates. In addition to numerical convergence 
(Sec. IIIB) and the neglect of the hyperfine interaction (Appendix B), we consider two 
additional sources of uncertainty, of which the first arises from neglecting the SE interaction 



[setting C = in Eq. (13)] and the second from inaccuracies in the ab initio interaction 
potential for the "^S^ state of N2 calculated in Sec. IIIA. 

To examine the sensitivity of the calculated dipolar relaxation rates to the SE interaction, 



we calculated the rates (A3) as functions of the SE parameter C (13). To estimate the range 
of variation of C, we used the ab initio results available for the ^SJ" electronic state. Figure 
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shows the PEC for the ^E+ state of N2 calculated for selected values of C using Eqs. (12) 



and (13). In the absence of the SE interaction, the ^S^ potential is identical to the ^E+ 
potential. As C increases, the ^S+ potential becomes deeper and shifts towards smaller 



R. While the Heisenberg exchange Hamiltonian (10) cannot accurately describe the shape 



of the ab initio PEC for the ^SJ" state [56], the long-range part of the curve is fairly well 



reproduced at C = 0.3 E^- Based on the comparison presented in Fig. 12, we choose to vary 
C in the range from to 0.5 E^ with a grid step of 0.05 E^. 



As shown in Fig. 13, the rate constants for dipolar relaxation (A3) do not vary strongly 
with C, except at the lowest temperature studied (0.1 K). The dependence of the calculated 
rates on C is not monotonous, and the largest deviation from C = values used as a reference 
in this work does not exceed 15 %. A similar lack of sensitivity has recently been observed in 
quantum calculations of dipolar relaxation in collisions of spin-polarized Eu atoms [7D] . We 
emphasize, however, that the SE interaction in the Eu2 collision complex is several orders 
of magnitude weaker than in N2, so the range of SE splittings probed in Ref. [70] was much 



narrower than explored in this work. The results presented in Fig. [13] therefore suggest that 
the rate constants for dipolar relaxation are insensitive to much larger variations of the SE 
interaction (on the order of several eV). 
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To understand how the calculated inelastic rates are affected by the uncertainties in the 
'''S+ interaction potential, we scaled the potential by a constant factor A and calculated 
the temperature dependence of (-ft"in) for 20 equally spaced values of A in the interval 0.9 
- 1.1 (A = 1 corresponds to unsealed potential B). The relatively small range of A chosen 
reflects the high level of accuracy of the ab initio interaction potential presented in Sec. IIIA 



(<10%). Figure 14 shows the dependence of (-ft'in(A)) for three selected temperatures. At 
T = 0.6 K, the variation of (K^^) is within 15 % for the whole range of A, demonstrating 
that our results are robust against both lessening (A < 1) and deepening (A > 1) of the 
interaction potential. The sensitivity of the calculated relaxation rates to A increases at 
low temperatures, reaching a maximum at T = 0.1 K. This is an expected result since the 
variation of the cross section with A should be most pronounced in the ultracold limit, where 
the s-wave scattering cross section exhibits a resonance-like variation as a function of A [71] . 



Figure 15 shows the temperature dependence of the calculated dipolar relaxation rates 
for ^^N + ^^N collisions. The error bars represent maximum possible deviations from the 
mean value of (-fCm) deflned as the value calculated in the absence of the SE interaction 
(C = 0) for unsealed potential B (A = 1). We evaluate the error bars by flnding the extrema 
of the calculated functions (-ft'in(A)) and (i^in(C)) for each T. The rate constants decrease 
and the error bars shrink with increasing temperature. The results presented in Table III 
indicate that while imperfections in the ^S+ interaction potential are the dominant source 
of uncertainty at temperatures below 0.6 K, omission of the SE interaction introduces the 
same amount of error at T = 0.6 K and becomes the major source of uncertainty above 



this temperature. From Figs. [13] and [14] we observe that scaling the interaction potential 
tends to increase the inelastic rates, whereas varying the strength of the SE interaction does 
not always lead to the monotonous variation. As a result, uncertainties in the interaction 
potential determine the upper error bar at T < 0.6 K, and those in the SE interaction 
determine the lower error bar at all temperatures. 

Figure 16 shows the variation of the elastic collision rates for ^^N + ^"^N with A. Strong 



sensitivity to A is apparent over the whole temperature range. As discussed in Sec. IIC, 
the calculated ratio of the rate constants for elastic scattering and dipolar relaxation in 
i4]\j _|_ i4]\j collisions is consistent with the measured value of 7 shown in Fig. [s] 

We note that the uncertainties arising from lack of knowledge of the SE interaction can 
be reduced by performing calculations with the accurate ab initio interaction potential for 
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the ^S+ electronic state of N2, that is, using the exact Eq. (10) instead of the approximate 



Heisenberg parametrization (11). Such calculations are currently in progress [72], and pre- 
liminary results show that the calculated relaxation rates are very close (within 5%) to the 
results obtained with the SE interaction omitted, thereby lending support to our claim that 
the major source of uncertainty in the calculated dipolar relaxation rates comes from the 
interaction potential for the ^S+ electronic state, rather than the SE interaction. 

IV. SUMMARY 

We have presented a combined experimental and theoretical study of collisional properties 
of cold spin-polarized atomic nitrogen. We have trapped large numbers of ^'^N and ^^N atoms 
for tens of seconds and measured their dipolar relaxation rates at 600 mK. Based on these 
measurements and theoretical calculations of trap loss rates, we have determined the number 
density of trapped N atoms to be (5±2) x 10^^ cm~^. 

Our theoretical analysis of dipolar relaxation in N + N collisions is based on accurate ab 
initio interaction potentials for the ^E+ state of N2 computed using highly correlated coupled 
cluster methods (Sec. IIIA). By solving the multichannel scattering problem, we obtained 
the cross sections and rate constants for dipolar relaxation in N + N collisions over a wide 
range of collision energies and magnetic fields (Sec. IIIB). The calculated relaxation rates 
for both ^^N and ^^N isotopes are similar in the multiple partial wave regime (T > 5 mK), 
but display a very different behavior at ultralow temperatures due to the effects of quantum 



statistics (Figs. ^ and 10). The rate constants for dipolar relaxation in N + N collisions 
are on the order of 10^^"^ cm^/s, indicating that spin-polarized N atoms are stable against 
collisional relaxation in the temperature range between 1 mK and 1 K. The results presented 
in Fig. [To] indicate that sympathetic cooling of paramagnetic molecules with N atoms will 
be efficient provided the probabilities for inelastic relaxation in N-molecule collisions are not 
very large. At T < 1 mK, the elastic cross section for ^^N + ^^N decreases dramatically 



and the elastic-to-inelastic ratio for ^^N + ^^N drops below 100 (Fig. 10). Thus, neither 
^^N nor ^^N appears suitable for sympathetic cooling of molecules below 1 mK. It might be 
possible to further reduce the temperature of trapped molecules via evaporative cooling at 
low magnetic fields [67] once N atoms are removed from the trap. 

In agreement with a recent theoretical study of dipolar relaxation in Eu + Eu collisions 
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[70j . we found that the calculated rate constants for trap loss in collisions of spin-polarized N 
atoms are not sensitive to the magnitude of the SE interaction. We identified inaccuracies in 
the interaction potential for the ^S^ electronic state of N2 as the major source of uncertainty 
in our theoretical results. These inaccuracies (on the order of 10%) lead to large variations 



of the calculated relaxation rates at temperatures below 0.1 K (Fig. 14), but have a minor 



effect at the experimental temperature of 0.6 K (Fig. 15 and Table III), enabling accurate 
calibration of the trapped N atom density (Sec. II). The results shown in Fig. 15 demonstrate 
that rigorous quantum scattering calculations based on ab initio interaction potentials are 
capable of providing quantitative accuracy required for the interpretation of cold collision 
experiments in the temperature range between 0.1 and 0.7 K. 



The calculated rate constants for dipolar relaxation at temperatures below 0.1 K (Fig. 10) 



are subject to large uncertainties arising from imperfections in ab initio interaction poten- 



tials. As shown in Fig. [T5| the calculated trap loss rates at T = 0.1 K are only accurate to 
within a factor of 3. At lower temperatures, scattering resonances similar to those shown in 
Fig. [8] may have a profound effect on collision dynamics. The positions and widths of these 
resonances are extremely sensitive to tiny variations in the interaction potentials and hence 
cannot be predicted quantitatively. As in the case of ultracold collisions of alkali-metal atoms 
[llj . empirical adjustment of the ab initio interaction potentials may be required to obtain 
quantitative agreement with future experimental studies of N + N collisions at temperatures 
below 0.1 K. 

Our findings indicate that spin-polarized nitrogen atoms may have favorable coUisional 
properties over a wide range of temperatures and magnetic fields, making them promising 
candidates for future experiments on sympathetic cooling of open-shell molecules such as NH 
[T7] to temperatures ~1 mK. A detailed study of cold N -|- NH collisions in a magnetic trap 
will be presented in future work [73]. The moderate magnetic moment of N atoms (3/is) 
is large enough to enable efficient magnetic trapping and evaporative cooling [T71 |2S] and 
small enough to make collision-induced dipolar relaxation inefficient. The latter property is 
particularly important since large inelastic loss rates recently observed in collisions of highly 
magnetic atoms [TH [75j make these atoms unsuitable for sympathetic cooling of molecules 
in permanent magnetic traps. 
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Appendix A: Rate equations 

Here, we present the derivation of the rate equation dsl). Assuming that the process of 
trap loss is irreversible, the time decay of trapped N atom density n due to dipolar relaxation 
in binary N -|- N collisions can be described by the following rate equation pTl [65] 

-n= Y. ^M's^M'^^K^l {B,T)n\ (Al) 



where Kzz , , {B, T) are state-resolved rate constants for dipolar relaxation (21 ). The 

2 2^^Sa*^Sb 1 I 

weighting factors in Eq. ( |A1[ ) serve to distinguish between single spin-flip collisions, in which 
only one atom is lost {wi3 = 1) and double spin-flip collisions, in which both atoms are lost 
{wii = 2). Taking into account only the dominant relaxation channels shown in Fig. 
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we can rewrite the rate equation (Al) in the form 

-n = KUB,T)n\ (A2) 

where 



K,^iB, T) = \ [2K^ 3 , 1 1 (5, T) + K3 3 , 1 3 (fi, T)] (A3) 

22^22 22^22 

is the total rate constant for trap loss and the factor of 1/2 is introduced to account for indis- 



tinguishability of collision partners |65j- The right-hand side of Eq. (A3) can be evaluated 



in terms of the partial rate constants given by Eq. (21 ). 



The rate constants defined by Eq. (A3) characterize the dynamics of dipolar relaxation 



in the presence of a uniform magnetic field. The trapping field generated in our apparatus 
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(Fig. 1) is, however, highly inhomogeneous, so the calculated loss rates (A3) should be 
averaged over the magnetic field distribution of the trap. To do this, we assume a trapped 
sample density distribution of the form 

n{r) = noU{r,T) (A4) 

where f/(r, T) = exp[— fiB (r) / ksT] is the magnetic field distribution of the trap, no is 
the density of N atoms at the trap center, B{r) is the trapping field, and T is the atom 
temperature. The trapping field is calculated numerically from the known electromagnetic 
coil profiles, and then fit to an 11 term polynomial 



Integration of Eq. (A2) over the trap volume using the density distribution of Eq. (A4) 



where 



yields an expression for total trap loss 

is the average rate constant for trap loss, and the value of 1/7.6 comes from the numeric 
evaluation of the expression r J^^^ j-)^^ ~ yg for the experimental trap geometry (Fig. 1) [26] . 



Appendix B: Hyperfine interaction 

In order to justify the approximation of neglecting the hyperfine structure we made in 
Sec. IIIA, we performed test calculations of ^^N + ^^N collisions with the hyperfine structure 
included. The results for Ec = 0.6 K and C = 0.5 Eh are shown in Fig. [17) The inelastic 
cross sections for N atoms colliding in the uppermost Zeeman state / (Fig. [T]) are identical 
to those calculated with the hyperfine structure omitted, as expected for the fully spin- 
polarized Zeeman states. When the atoms collide in partially polarized Zeeman states j 
or k at B < 10 mT, they can exchange spin angular momentum via the SE interaction 



(10). The cross sections for collision-induced SE relaxation are typically much larger than 



those for dipolar relaxation \ll\, so the inelastic cross sections increase by a factor of 50-100 



as shown in Fig. [T7j As B increases, the states j, k, and i converge to the same limit 
Ms = 3/2, and the inelastic cross sections decrease monotonically, approaching the same 
limiting value calculated without taking into account the hyperfine structure (Sec. IIIA). 
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At the temperature and trap depths for the experiments described in Sec. II, N atoms at 
fields below 10 mT account for less than 10"*^ of the total number of trapped atoms, and 
therefore do not make a significant contribution to the total trap loss rate. 



[1] L. D. Carr, D. DeMille, R. V. Krems, and J. Ye, New J. Phys. 11, 055049 (2009). 

[2] P. Rabl, D. DeMille, J. M. Doyle, M. D. Lukin, R. J. Schoelkopf, and P. Zoller, 

Phys. Rev. Lett. 97, 033003 (2006). 
[3] I. Block, Nature (London) 453, 1016 (2008). 
[4] R. V. Krems, Phys. Chem. Chem. Phys. 10, 4079 (2008). 
[5] B. Friedrich and J. M. Doyle, ChemPhysChem 10, 604 (2009). 

[6] J. M. Doyle, B. Friedrich, J. Kim and D. Patterson, Phys. Rev. A, 52, R2515 (1995). 
[7] E. Narevicius, A. Libson, C. G. Parthey, I. Chavez, J. Narevicius, U. Even, and M. G. Raizen, 

Phys. Rev. A 77, 051401 (R) (2008). 
[8] S. Y. T. van de Meerakker and G. Meijer, Faraday Discuss. 142, 113 (2009). 
[9] M. Schneh and G. Meijer, Angew. Chem. Int. Ed. 48, 6010 (2009). 
[10] L. D. van Buuren, C. Sommer, M. Motsch, S. Pohle, M. Schenk, J. Bayerl, P. W. H. Pinkse, 
and G. Rempe, Phys. Rev. Lett. 102, 033001 (2009); M. Motsch, L. D. van Buuren, C. 
Sommer, M. Zeppenfeld, G. Rempe, and P. W. H. Pinkse, Phys. Rev. A 79, 013405 (2009). 
[11] J. Weiner, V. S. Bagnato, S. Zilio, P. S. Julienne, Rev. Mod. Phys. 71, 1 (1999). 
[12] K. M. Jones, E. Tiesinga, P. D. Lett, and P. S. Julienne, Rev. Mod. Phys. 78, 483 (2006). 
[13] G. Modugno, G. Ferrari, G. Roati, R. J. Brecha, A. Simoni, and M. Inguscio, Science 294, 

1320 (2001). 
[14] S. Aubin, S. Myrskog, M. H. T. Extavour, L. J. Leblanc, D. Mckay, A. Stummer, and J. H. 

Thywissen, Nat. Phys. 2, 384 (2006). 
[15] P. G. Mickelson, Y. N. Martinez de Escobar, M. Yan, B. J. DeSalvo, and T. C. Kihian, 

Phys. Rev. A 81, 051601(R) (2010). 
[16] S. C. Doret, C. B. Connolly, W. Ketterle, and J. M. Doyle, Phys. Rev. Lett. 103, 103005 

(2009). 
[17] M. T. Hummon, W. C. CampbeU, H.-I Lu, E. Tsikata, Y. Wang, and J. M. Doyle, 
Phys. Rev. A 78, 050702(R) (2008). 

25 



[18] M. Lara, J. L. Bohn, D. Potter, P. Soldan, and J. M. Hutson, Phys. Rev. Lett. 97, 183201 

(2006); Phys. Rev. A 75, 012704 (2007). 
[19] M. Tacconi, E. Bodo, and F. A. Gianturco, Theor. Chem. Ace. 117, 649 (2007). 
[20] P. S. Zuchowski and J. M. Hutson, Phys. Rev. A 79, 062708 (2009). 
[21] T. V. Tscherbul and R. V. Krems, Phys. Rev. Lett. 97, 083201 (2006); T. V. Tscherbul, G. 

C. Groenenboom, R. V. Krems and A. Dalgarno, Faraday Discuss. 142, 127 (2009). 
[22] J. van Veldhoven, H. L. Bethlem, and G. Meijer, Phys. Rev. Lett. 94, 083001 (2005); M. 

Schneh, P. Lutzow, J. van Veldhoven, H. L. Bethlem, J. Kupper, B. Friedrich, M. Schleier- 

Smith, H. Haak, and G. Meijer, J. Phys. Chem. A 111, 7411 (2007). 
[23] D. DeMille, D. R. Glenn, and J. Petricka, Eur. Phys. J. D 31, 375 (2004). 
[24] A. O. G. Wallis and J. M. Hutson, Phys. Rev. Lett. 103, 183201 (2009). 
[25] P. Barletta, J. Tennyson, and P. F. Barker, New J. Phys. 11, 055029 (2009). 
[26] M. T. Hunimon, PhD thesis. Harvard University (2010). 
[27] S. Hensler, J. Werner, A. Griesmaier, P. O. Schmidt, A Gorlitz, T. Pfau, S. Giovanazzi, and 

K. Rzazewski, Appl. Phys. B 77, 765 (2003). 
[28] E. Tsikata, W. C. CampbeU, M. T. Hummon, H.-I Lu, and J. M. Doyle, New J. Phys. 12, 

065028 (2010). 
[29] J. G. E. Harris, and R. A. Michniak, S. V. Nguyen, N. Brahms, W. Ketterle, and J. M. Doyle, 

Europhys. Lett. 67, 198 (2004). 
[30] S. F. Adams and T. A. Miller, Chem. Phys. Lett. 295, 305 (1998). 

[31] G. J. Bengtsson, J. Larsson, S. Svanberg, and D. D. Wang, Phys. Rev. A 45, 2712 (1992). 
[32] D. J. Bamford, L. E. Jusinski, and W. K. Bischel, Phys. Rev. A 34, 185 (1986). 
[33] D. J. Bamford, A. P. Hickman, M. J. Dyer, and W. K. Bischel, J. Opt. Soc. Am. B. 5, 1369 

(1988). 
[34] M. G. Littman and H. J. Metcalf, App. Opt. 17, 2224 (1978). 
[35] T. T. Kajava, H. M. Lauranto, and R. R. E. Salomaa, App. Opt. 31, 6987 (1992). 
[36] L. A. Westling, M. G. Raymer, and J. J. Snyder, J. Opt. Soc. Am. B. 1, 150 (1984). 
[37] Moore, Charlotte E. Atomic energy levels as derived from the analyses of optical spectra. 

Volume 1, U.S. National Bureau of Standards, 1971. 
[38] J. M. Hirsch, G. H. Zimmerman, D. J. Larson, and N. F. Ramsey, Phys. Rev. A 16, 484 

(1977). 

26 



[39] W. Ketterle and N. J. VanDruten, Adv. At. Mol. Opt. Phys. 37, 181 (1996); R. deCarvalho 

and J. M. Doyle, Phys. Rev. A 70, 053409 (2004). 
[40] K. Omidvar, Phys. Rev. A 22, 1576 (1980). 
[41] T. H. Dunning, Jr., J. Chem. Phys. 90, 1007 (1989) 

[42] R. A. Kendall, T. H. Dunning, Jr. and R. J. Harrison, J. Chem. Phys. 96, 6796 (1992) 
[43] P. J. Knowles, C. Hampel, and H.-J. Werner, J. Chem. Phys. 99, 5219 (1993); Erratum: J. 

Chem. Phys. 112, E3106-3107 (2000). 
[44] MOLPRO, version 2008.1, a package of ab initio programs, H.-J. Werner, P. J. Knowles, R. 



Lindh, F. R. Manby, M. Schiitz, and others, see http://www.molpro.net 
[45] D. Feller and J. A. Sordo, J. Chem. Phys. 112, 5604 (2000). 
[46] T. S. Ho and H. Rabitz, J. Chem. Phys. 104, 2584 (1996). 



[47] See supplementary material at |http://link. aps.org/supplemental/10.1103/PhysRevA.82.042718 



for the Fortran codes used in this work to generate potentials A and B as functions of R. 
[48] K. A. Peterson, D. E. Woon, and T. H. Dunning, Jr., J. Chem. Phys. 100, 7410 (1994). 
[49] MRCC, a string-based quantum chemical program suite written by M. Kallay. For more in- 



formation see Ref. [50] and http://www.mrcc.hu 



[50] M. Kallay and P. R. Surjan, J. Chem. Phys. 115, 2945 (2001); M. Kallay, P. G. Szalay and P. 

R. Surjan, J. Chem. Phys. 117, 980 (2002); M. Kallay, J. Gauss and P. G. Szalay, J. Chem. 

Phys. 119, 2991 (2003); M. Kallay and J. Gauss, J. Chem. Phys. 120, 6841 (2004); M. Kallay 

and J. Gauss, J. Chem. Phys. 121, 9257 (2004); M. Kallay and J. Gauss, J. Chem. Phys. 123, 

214105 (2005). 
[51] H. Partridge, S. R. Langhoff, and C. W. Bauschlicher, Jr., J. Chem. Phys. 84, 6901 (1986). 
[52] X. Chu and A. Dalgarno, J. Chem. Phys. 121, 4083 (2004). 
[53] L. Anderson, F. M. Pipkin, and J. C. Baird, Phys. Rev. 116, 87 (1959). 
[54] H. T. C. Stoof, J. M. V. A. Koelman, and B. J. Verhaar, Phys. Rev. B 38, 4688 (1988). 
[55] X. Li and J. Paldus, J. Chem. Phys. 129, 054104 (2008); X. N. Zhang, D. H. Shi, J. F. Sun, 

and Z. L. Zhu, Chin. Phys. B 19, 013501 (2010). 
[56] H. Partridge, S. R. Langhoff, C. W. Bauschlicher, Jr., and D. W. Schwenke, J. Chem. Phys. 88, 

3174 (1988). 
[57] W. T. Zemke and W. C. Stwalley, J. Chem. Phys. Ill, 4962 (1999). 



27 



[58] R. Cote, M. Jamieson, Z-C. Yan, N. Geum, G.-H. Jeung, and A. Dalgarno, Phys. Rev. Lett. 84, 
2806 (2000). 

[59] B. M. Smirnov and M. I. Chibisov, Sov. Phys. JETP 21, 624 (1965). 

[60] R. N. Zare, Angular momentum (Wiley, New York, 1988). 

[61] Z. Li and R. V. Krems, Phys. Rev. A 75, 032709 (2007). 

[62] Heavy atoms in fully spin-polarized Zeeman states such as ^^^Cs(-F = 3, mp = 3) undergo 
inelastic relaxation induced by a second-order spin-orbit interaction, see P. J. Leo, E. Tiesinga, 
P. S. Julienne, D. K. Walter, S. Kadlecek, and T. G. Walker, Phys. Rev. Lett. 81, 1389 (1998). 

[63] B. Zygelman and A. Dalgarno, J. Phys. B 35, L441 (2002); B. Zygelman, 



|arXiv:physics/0209052| (2002). 
[64] J. M. Hutson, E. Tiesinga, and P. S. Juhenne, Phys. Rev. A 78, 052703 (2008). 
[65] J. P. Burke, Jr., PhD thesis. University of Colorado (1999). 
[66] R. V. Krems and A. Dalgarno, J. Chem. Phys. 120, 2296 (2004). 
[67] T. V. Tscherbul, Yu. V. Suleimanov, V. Aquilanti, and R. V. Krems, New J. Phys. 11, 055021 

(2009). 
[68] D. E. Manolopoulos, J. Chem. Phys. 85, 6425 (1986). 

[69] Z. Pavlovic, R. V. Krems, R. Cote, and H. R. Sadeghpour, Phys. Rev. A 71, 061402(R) (2005). 
[70] Y. V. Suleimanov, Phys. Rev. A 81, 022701 (2010). 
[71] M. Cvitas, P. Soldan, J. M. Hutson, P. Honvault, and J. M. Launay, J. Chem. Phys. 127, 

074302 (2007). 
[72] J. Klos, to be published. 
[73] M. T. Hummon, T. V. Tscherbul, J. Klos, H-I Lu, E. Tsikata, W. C. Campbell, A. Dalgarno, 

and J. M. Doyle, Phys. Rev. Lett. 106, 053201 (2011). 
[74] S. V. Nguyen, R. deCarvalho, and J. M. Doyle, Phys. Rev. A 75, 062706 (2007). 
[75] C. B. Connolly, Y. S. Au, S. C. Doret, W. Ketterle, and J. M. Doyle, Phys. Rev. A 81, 

010702(R) (2010). 



28 



Table I. Parameters and typical values for atomic nitrogen excitation. 



description 



symbol 



typical value units 



two-photon cross section [iO] 

excitation pulse energy 

beam waist 

effective excitation length 

pulse duration (FWHM) 

resonant line-shape value g{0) (2/7r)(27r x 10 GHz)"^ s 

2nd-order photon correlation coefficient G^^\0) 2 



E 

Wo 

'ex 



1.37 10-36 cm^ 
^ 0.6 mJ 

120 jjm 
2 mm 

9.5 ns 



photon collection efficiency 



a 



10" 
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Table II. Ro-vibrational levels of N2('^S+) supported by potentials A and B. The level 
energies are given in cm"-"^ relative to the N(^S'3/2) + N(*^S'3/2) dissociation limit in the 
absence of a magnetic field. The magnetic dipole interaction is not included in the bound- 
state calculations. 



v,(i 


i^Nj 




^^N2 






Potential A 


Potential B 


Potential A 


Potential B 


0, 


-17.38 


-19.10 


-17.72 


-19.47 


0,1 


-17.08 


-18.80 


-17.45 


-19.19 


0,2 


-16.50 


-18.21 


-16.90 


-18.63 


0, 3 


-15.62 


-17.32 


-16.08 


-17.80 


0,4 


-14.46 


-16.14 


-14.98 


-16.69 


0, 5 


-13.01 


-14.66 


-13.62 


-15.30 


0, 6 


-11.28 


-12.91 


-12.00 


-13.66 


0,7 


-9.28 


-10.88 


-10.12 


-11.74 


0, 8 


-7.02 


-8.57 


-7.99 


-9.58 


0, 9 


-4.51 


-6.01 


-5.62 


-7.16 


0, 10 


-1.76 


-3.20 


-3.03 


-4.51 


0, 11 


- 


-0.162 


-0.221 


-1.64 


1,0 


-3.61 


-4.30 


-4.02 


-4.75 


1,1 


-3.41 


-4.09 


-3.83 


-4.56 


1,2 


-3.01 


-3.68 


-3.45 


-4.16 


1,3 


-2.43 


-3.08 


-2.87 


-3.58 


1,4 


-1.67 


-2.29 


-2.15 


-2.82 


1,5 


-0.75 


-1.32 


-1.25 


-1.88 


1,6 


- 


-0.205 


-0.221 


-0.79 


2,0 


-0.066 


-0.125 


-0.129 


-0.217 


2,1 


-0.011 


-0.055 


-0.062 


-0.137 
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Table III. The calculated rate constants (-ft'in) for ^^N (in units of 10~^^ cm^/s) vs temper- 
ature (in K). The maximum relative uncertainties with respect to the mean value calculated 
for potential B, C = and A = 1 are presented in the third column. The error bars are 
based on two sets of calculations using (i) A = 0.90, . . . , 1.10 with a grid spacing of 0.01 for 
C = and (ii) C = 0, . . . , 0.5 Eh with a grid spacing of 0.05 Eh for A = 1. Also indicated 
are the dominant sources of uncertainty in the calculated rates arising from inaccuracies in 
the '^S+ interaction potential (IP) and omission of the SE interaction (SE). 



Temperature 


(i^in) 


Uncertainty (%) 


Source 


0.1 


3-8-0.6 


184.2 


IP 


0.2 


A c,+2.6 


57.8 


IP 


0.3 


4 3+^-^ 

^-<J-0.5 


25.6 


IP 


0.4 


4-0_o.5 


17.5 


IP 


0.5 


q 7+0.5 
-J-' -0.4 


13.5 


IP 


0.6 


q /I+0.4 


11.8 


IP and SE 


0.7 


q 9+0.2 


12.5 


SE 
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Fig. 1. Diagram of trapping apparatus. 
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Fig. 2. Trapped nitrogen spectrum, fitted to a Lorentzian profile. 
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Fig. 3. Time decay of trapped nitrogen. Initial trap loading occurs at t = 0. The trap loss is fitted 
to the solution of Eq. [S] to determine the trap loss parameters. 
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Fig. 4. (i^ioss) vs r/. The solid line is a fit of the data to Eq. 5 The dashed line is a fit of the data 
to [5] with 7 set to the value calculated in Sec. Ill, 7theory = 1000. 
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Fig. 5. The calculated (crosses) and measured (circles) ratios 7 = (i^ci)/(-f^'in) for -"^"^N. The theo- 
retical ratios and error bars are calculated based on the data shown in Figs. 15 and 16 as described 
in Sec. IIIC. 
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Fig. 6. Ab initio interaction potentials for the S^ electronic state of N2 calculated as functions 
of the internuclear distance. Dashed line - potential A computed using the UCCSD[T] method, 
full line - potential B computed using the MRCCSDT method. 
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Fig. 7. Zeeman energy levels of ^^N as functions of the applied magnetic field. The inset shows 
details of the hyperfine structure at low fields. The states are labeled according to their spin 
projections Msp^- The highest-energy low-field-seeking state \Sp^ = 2>/2,Msj^ = 3/2) is shown by 
the dashed line. 
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Fig. 8. Cross sections for dipolar relaxation (full lines) and elastic scattering (dashed lines) in ^^N 
+ -"^^N and -"^^N + -"^^N collisions plotted vs collision energy Ec for i? = 0.1 T and C = 0. The peak 
around Ec = 1.3 cm^^ corresponds to an £ = 7 shape resonance. 
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Fig. 9. (a) Cross sections for dipolar relaxation in ^^N + ^^N collisions as functions of collision 
energy (in cni^^) and magnetic field (in T). Note the presence of the 1 = 7 shape resonance marked 
in Fig. 8 and its evolution with magnetic field, (b) Same but for -"^^N + ^^N collisions. The cross 
sections are evaluated for C = 0. 
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Fig. 10. (a) Rate constants for dipolar relaxation in collisions of ^^N atoms (circles) and ^^N atoms 
(triangles) as functions of temperature, (b) Thermally averaged ratios of the rate constants for 
elastic scattering and dipolar relaxation. The magnetic field is 0.1 T. 



41 



1x10 



T = 0.6K 



5x10 



Total 



3 

a 

X 






8x10 



4x10 




0.5 1 1.5 

Magnetic field (T) 



Fig. 11. (a) Partial rate constants for dipolar relaxation in ^^N + ^^N collisions as functions of 
magnetic field calculated with potential A for (a) T = 0.6 K, (b) T = 0.1 K and C = 0. The total 
rate constant for dipolar relaxation is shown by tlie dashed line. 
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Fig. 12. Potential energy curve for the j4'^S^ electronic state of N2 calculated for different values 
of the SE parameter C indicated in the graph. The MRCI+Q results of Partridge et al. [IS] are 
shown by triangles. 
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Fig. 13. Rate constants for dipolar relaxation (A6) calculated for N + N as functions of the 



SE parameter C for T = 0.6 K and C = 0. A = 1 corresponds to unsealed potential B. 



44 



10 



§ 6 









1 
Temperature (K) 


— 


- \ 


^ 




•-•0.1 
♦-♦0.3 
^—^0.6 


/ 


- 




^ — • — • — 1 


1 











0.9 



0.95 1 1.05 

Potential scaling parameter ^ 



1.1 



Fig. 14. Rate constants for dipolar relaxation (A6) calculated for N + N as functions of the 



potential scaling parameter A for T = 0.6 K and C = 0. A = 1 corresponds to unsealed potential B. 
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Fig. 15. Rate constants for dipolar relaxation (A6) calculated for N + N as functions of 



temperature. The error bars are calculated as explained in the text (see also Table III). 
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Fig. 16. Elastic collision rates for ^ N + ^^N averaged over the magnetic field distribution of the 
trap as functions of temperature. 
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Fig 17. Cross sections for inelastic relaxation in collisions of ^^N atoms in the low-field-seeking 
hyperfine states j'j, kk, and // (full lines) and jk, kl, jl (dashed lines) calculated for Ec = 0.6 K and 
C = 0.5 Eh- Also shown are the cross sections calculated with the hyperfine interaction omitted 
(triangles) for the \Msj^ = 3/2)\Ms^ = 3/2) initial channel. The hyperfine levels are labeled as 
shown in the inset of Fig. [7| 
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